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EXPLOITATION  OF  MULTIPLE  SOLUTIONS  OF  THE 
NAVIER-STOKES  EQUATIONS  TO  ACHIEVE  RADICALLY  IMPROVED  FLIGHT 

Objectives  of  the  Research 

It  is  generally  accepted  that  all  fluid  flow  phenomena  of  practical  use, 
including  turbulence,  can  be  represented  by  the  Navier-Stokes  equations.  The 
Navier-Stokes  equations  are  a  set  of  coupled  nonlinear  partial  differential  equations 
representing  the  conservation  of  mass,  momentum,  and  energy.  The  nonlinearity  of 
the  equations  allows  the  representation  of  turbulence,  shock  waves,  and  slip  surfaces. 
If  the  equations  are  Reynolds  averaged,  then  the  ability  to  represent  the  details  of 
turbulence  is  lost,  but  the  nonlinearity  necessary  to  represent  shock  waves  and  slip 
surfaces  is  retained.  Little  is  known  about  the  existence  and  uniqueness  of  the 
solutions  to  nonlinear  equations.  A  common  test  of  existence  is  to  see  if  the 
solutions  of  the  Navier-Stokes  equations  closely  model  the  behavior  observed  in  wind 
tunnel  tests.  It  is  reasonable  to  speculate  that  if  the  Navier-Stokes  equations  had 
several  solutions,  one  of  which  closely  models  an  observed  flow,  then  the  other 
solutions  must  be  physically  realizable,  although  probably  unstable.  This  raises  the 
intriguing  question  of  whether  these  realizable  but  unstable  solutions  represent  a 
radical  change  in  flow  behavior  that  may  be  of  great  benefit  in  a  practical 
aerodynamic  sense. 

In  the  present  work,  the  conventional  solution  is  defined  as  the  solution  that 
would  be  expected  from  an  extrapolation  from  previous  knowledge  and  a  "phantom'' 
solution  is  defined  as  any  additional  solution. 
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One  type  of  phantom  solution  that  is  found  in  potential  transonic  flow  theory, 
a  simpler  subset  of  the  Navier-Stokes  equations,  indicates  that  there  are  several 
different  shock  patterns  that  satisfy  the  boundary  conditions.  The  typical  example 
of  this  type  of  phantom  solution  is  that  of  the  flow  around  a  symmetric  airfoil  at 
zero  angle  of  attack  giving  lift;  the  phenomena  appears  also  for  conventional  lifting 
airfoils.  An  example^  of  the  variation  of  lift  coefficient,  Cl,  with  angle  of  attack 
a,  for  a  RAE2822  airfoil  is  shown  in  Figure  1.  It  may  be  observed  that  at  a  = 
0.6  degrees,  Cl  can  have  a  value  of  about  0.5  or  about  1.4.  Phantom  solutions 
have  also  been  found  for  three-dimensional  flows. ^  This  phenomena  was  first 
reported  by  Steinhoff  and  Jameson.^  These  phantom  solutions  appear  when  some 
degree  of  asymmetry,  either  in  the  algorithm  or  in  the  initial  conditions,  is 
introduced.  The  asymmetry  may  also  be  introduced  into  the  solution  process  by  a 
physical  angle  of  attack.  It  appears  that  there  is  a  certain  range  of  Mach  number 
in  which  the  phantom  solution  is  preferred  to  the  conventional  solution.  These 
phantom  solutions  have,  until  recently,  not  been  observed  in  equation  sets  other 
than  the  potential  equation;  although,  it  has  been  speculated  by  Williams  et  al^ 
and  Nixon^  that  they  would  exist  for  the  Euler  equations.  Nixon^  has  recently 


computed  phantom  solutions  for  an  equation,  the  "TSD-Euler"  equation,  which 
closely  resembles  the  Euler  equations.  For  configurations  that  would  normally  give 
lift,  the  presence  of  one  of  the  phantom  solutions  can  increase  the  lift  considerably, 
by  a  factor  of  two  or  three  (see  Figure  1).  It  is  therefore  a  more  than  academic 
point  to  consider  whether  these  phantom  solutions  can  exist  in  practice. 
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The  present  work  is  concerned  with  a  continuing  investigation  of  the 
appearance  of  phantom  solutions  for  the  TSD-Euler  equation  with  a  view  to 
determining  if  such  solutions  can  exist  for  the  Navier-Stokes  equations  and  thus 
indicating  that  they  may  be  physically  realizable.  The  evidence  gained  during  the 
present  study  does  indicate  that  phantom  solutions  could  occur  in  reality  and  thus 
may  introduce  a  new  aspect  of  aerodynamics,  namely  the  exploitation  of  nonlinear 
aerodynamic  phenomena. 

It.  was  anticipated  that  the  tasks  that  were  to  be  performed  during  the  course 
of  the  contract  should  be  as  follows: 

(a)  Study  further  the  occurrence  of  phantom  solutions  of  the  TSD-Euler  equation 
to  get  a  clearer  idea  of  the  mechanisms  and  control  of  such  solutions. 
Attempt  to  develop  a  means  of  differentiating  a  phantom  solution  from  a 
conventional  solution. 

(b)  Extend  the  XTRAN2L  code  with  boundary  layer  to  treat  the  TSD-Euler 
equation.  Determine  the  effect  of  a  boundary  layer  on  the  occurrence  and 
control  of  phantom  solutions. 

(c)  Using  the  ARC2D  code,  attempt  to  reproduce  the  phantom  solutions  for  the 
Euler  equations. 

(d)  Using  the  ARC2D  code,  attempt  to  reproduce  the  phantom  solutions  for  the 
Navier-Stokes  equations. 
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Status  of  Research  Effort 


The  TSD-Euler  equation,  derived  in  previous  work  was  examined  for  phantom 
solutions  in  an  attempt  to  determine  the  effect  on  entropy  and  vorticity  on  the 
appearance  of  such  solutions.  It  was  found  that  phantom  solutions  can  be  found 
for  this  equation  although  for  a  different  range  of  Mach  numbers  than  for  the  TSD 
equation,  at  least  for  the  airfoils  considered.  Following  earlier  work  a  source  of 
vorticity  was  introduced  on  the  airfoil  surface  and  this  enhanced  the  probability  of 
the  appearance  of  phantom  solutions.  Closer  study  indicated  that  the  dominant 
effect  of  the  vorticity  introduction  is  to  change  the  effective  geometry  of  the  airfoil 
with  the  implication  that  an  airfoil  could  be  designed  to  give  a  phantom  solution  to 
the  TSD-Euler  equation.  The  vorticity  introduction  can  be  regarded  as  a  method 
of  flow  control  which  mimics  the  effect  of  a  change  in  geometry. 

A  test  of  when  a  phantom  solution  can  appear  has  been  developed  and  it 
indicates  that  the  dominant  feature  determining  the  appearance  of  phantom  solutions 
is  the  local  Mach  number  distribution;  the  entropy  or  vorticity  in  the  flow  field 
plays  only  a  secondary  role.  The  numerical  verification  of  the  test  has  some 
problems  due  principally  to  the  numerical  inaccuracies  in  the  numerical  solutions  of 
the  flow  and  in  the  evaluation  of  the  test  parameters.  The  test  indicates  that  a 
phantom  solution  can  be  initiated  by  designing  the  airfoil  to  give  a  suitable 
distribution  of  Mach  number  in  the  flow  field. 
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The  inference  of  the  study  of  the  TSD-Euler  equation  is  that  entropy  and 
vorticity  do  not  inhibit  the  appearance  of  phantom  solutions  but  simply  change  the 
freestream  conditions  for  which  these  solutions  can  appear.  Since  vorticity  does  not 
inhibit  the  appearance  of  the  phantom  solutions  it  is  pertinent  to  ask  if  the  Navier- 
Stokes  equations  have  phantom  solutions.  It  is  known  that  for  a  small  Mach 
number  range  that  the  Navier-Stokes  equations  do  have  multiple  solutions,  are 
symmetric  and  are  oscillatory.  This  behavior  is  shown  to  be  explainable  b>  the 
idea  of  phantom  solutions,  although  this  may  not  be  the  only  explanation.  The 
multiple  solutions  of  the  Navier-Stokes  equations  are  apparent  in  wind  tunnel  tests, 
indicating  their  realism.  The  question  that  now  arises  is  whether  phantom  solutions 
can  be  made  to  exist  stably  in  a  real  environment  thus  introducing  a  new  area  of 
study  for  aerodynamicists  and,  possibly,  introducing  a  new  mechanism  for  flight.  In 
order  to  achieve  this  the  phantom  solutions  must  be  stable  and  some  form  of  flow 
control  is  necessary. 

If  the  test  derived  in  this  present  work  can  be  made  more  accurate 
numerically  it  should  be  possible  to  design  an  airfoil  to  give  phantom  solutions  over 
a  range  of  flight  conditions.  Again,  some  form  of  flow  control  is  necessary.  If 
such  an  airfoil  can  be  designed  then  phantom  solutions  to  the  Navier-Stokes 
equations  may  be  found,  together  with  some  means  of  stabilizing  such  solutions. 
After  this  stage  a  wind  tunnel  test  could  be  conducted.  To  achieve  such  an  end  it 
is  necessary  to  learn  much  more  about  the  nature  of  all  phantom  solutions,  not  just 
those  considered  here,  that  is,  the  transonic  solutions.  A  more  complete  discussion 
of  this  work  is  contained  in  the  Appendix. 
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It  proved  difficult  to  run  the  XTRAN2L  code  with  the  boundary  layer.  It 
was  decided  not  to  attempt  to  compute  Euler  solutions  because  of  the  unrealistic 
nature  of  the  boundary  conditions.  Phantom  solutions  to  the  Navier-Stokes 
equations  were  computed  but  are  unsteady.  Efforts  to  stabilize  these  solutions 
failed. 

(c)  Written  Publications 

Further  Study  of  Phantom  Solutions  to  the  TSD-Euler  Equation 
(to  be  submitted  for  publication  in  Acta  Mechanica) 

(d)  Personnel 

Dr.  David  Nixon 

Dr.  Steven  C.  Caruso 

Dr.  Mohammad  Farshchi 

No  degrees  have  been  awarded  as  a  result  in  this  work. 

(e)  Interactions 

The  work  has  been  discussed  with  Dr.  Joseph  Shang,  Dr.  James  Olsen,  and 
Dr.  Jan  Lee  at  AFWAL.  No  other  interactions  have  taken  place. 

(f)  No  new  discoveries  or  inventions,  patent  disclosures  or  specific  applications 
have  stemmed  from  the.  research  other  than  those  discoveries  detailed  in  this  report. 
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(g)  None. 
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APPENDIX 


FURTHER  STUDY  OF  PHANTOM  SOLUTIONS 
TO  THE  "TSD-EULER"  EQUATION 


Introduction 


It  is  generally  accepted  that  all  fluid  flow  phenomena  of  practical  use, 
including  turbulence,  can  be  represented  by  the  Navier-Stokes  equations.  The 
Navier-Stokes  equations  are  a  set  of  coupled  nonlinear  partial  differential  equations 
representing  the  conservation  of  mass,  momentum,  and  energy.  The  nonlinearity  of 
the  equations  allows  the  representation  of  turbulence,  shock  waves,  and  slip  surfaces. 
If  the  equations  are  Reynolds  averaged,  then  the  ability  to  represent  the  details  of 
turbulence  is  lost,  but  the  nonlinearity  necessary  to  represent  shock  waves  and  slip 
surfaces  is  retained.  Little  is  known  about  the  existence  and  uniqueness  of  the 
solutions  to  nonlinear  equations.  A  common  test  of  existence  is  to  see  if  the 
solutions  of  the  Navier-Stokes  equations  closely  models  the  behavior  observed  in 
wind  tunnel  tests.  It  is  reasonable  to  speculate  that  if  the  Navier-Stokes  equations 
had  several  solutions,  one  of  which  closely  models  an  observed  flow,  then  the  other 
solutions  must  be  physically  realizable,  although  probably  unstable.  This  raises  the 
intriguing  question  of  whether  these  realizable  but  unstable  solutions  represent  a 
radical  change  in  flow  behavior  that  may  be  of  great  benefit  in  a  practical 
aerodynamic  sense. 

In  this  paper  that  follows,  the  conventional  solution  is  defined  as  the  solution 
that  would  be  expected  from  an  extrapolation  from  previous  knowledge  and  a 
“phantom"  solution  is  defined  as  any  additional  solution. 
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One  type  of  phantom  solution  that  is  found  in  potential  transonic  flow  theory, 
a  simpler  subset  of  the  Navier-Stokes  equations,  indicates  that  there  are  several 
different  shock  patterns  that  satisfy  the  boundary  conditions.  The  typical  example 
of  this  type  of  phantom  solution  is  that  of  the  flow  around  a  symmetric  airfoil  at 
zero  angle  of  attack  giving  lift;  the  phenomena  appears  also  for  conventional  lifting 
airfoils.  An  example^  of  the  variation  of  lift  coefficient,  C^,  with  angle  of  attack 
a,  for  a  RAE2822  airfoil  is  shown  in  Figure  1.  It  may  be  observed  that  at  a  = 
0.6  degrees,  can  have  a  value  of  about  0.5  or  about  1.4.  Phantom  solutions 
have  also  been  found  for  three-dimensional  flows. ^  This  phenomena  was  first 
reported  by  Steinhoff  and  Jameson. ^  These  phantom  solutions  appear  when  some 
degree  of  asymmetry,  either  in  the  algorithm  or  in  the  initial  conditions,  is 
introduced.  The  asymmetry  may  also  be  introduced  into  the  solution  process  by  a 
physical  angle  of  attack.  It  appears  that  there  is  a  certain  range  of  Mach  number 
in  which  the  phantom  solution  is  preferred  to  the  conventional  solution.  These 
phantom  solutions  have,  until  recently,  not  been  observed  in  equation  sets  other 
than  the  potential  equation;  although,  it  has  been  speculated  by  Williams  et  al^ 
and  Nixon^  that  they  would  exist  for  the  Euler  equations.  Nixon^  has  recently 
computed  phantom  solutions  for  an  equation,  the  "TSD-Euler"  equation,  which 
closely  resembles  the  Euler  equations.  For  configurations  that  would  normally  give 
lift,  the  presence  of  one  of  the  phantom  solutions  can  increase  the  lift  considerably, 
by  a  factor  of  two  or  three  (see  Figure  l).  It  is  therefore  a  more  than  academic 
point  to  consider  whether  these  phantom  solutions  can  exist  in  practice. 
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The  present  paper  is  concerned  with  a  continuing  investigation  of  the 
appearance  of  phantom  solutions  for  the  TSD-Euler  equation  with  a  view  to 
determining  if  such  solutions  can  exist  for  the  Navier-Stokes  equations  and  thus 
indicating  that  they  may  be  physically  realizable.  The  evidence  gained  during  the 
present  study  does  indicate  that  phantom  solutions  could  occur  in  reality  and  thus 
may  introduce  a  new  aspect  of  aerodynamics,  namely  the  exploitation  of  nonlinear 
aerodynamic  phenomena. 


Basic  Equations 

The  TSD-Euler  equation  [5]  is  derived  from  the  classic  Transonic  Small  Disturbance 
(TSD)  equation  and  is  given  by 

(  A  -  k  Z  )  ♦  »y  -  (  f  )x  vi) 

x 

where  u  and  v  are  perturbation  velocities  in  the  x  and  y  directions,  respectively, 
nondimensionalized  with  respect  to  freestream  values.  S  is  the  entropy  change,  R  is 
the  gas  constant  and 


P2  =  1  -  k  =  (7  +  1)  M^'5  (2) 

where  M«  is  the  freestream  Mach  number.  If  ^  is  a  perturbation  velocity  potential, 
and  ^  is  a  vector  potential,  then 


A-3 


and 

4  +  4  =  -u  (4) 

rxx  ryy  v  ' 

where  the  vorticity  u  is  related  to  the  entropy  S  by  Crocco's  theorem.  The 
entropy  generated  ^  at  a  shock  is  given  by 

|  ~  - To-—  (M+2  _  l)3  (5) 

R  (7  +  l)2 

where  M+  is  the  Mach  number  just  ahead  of  the  shock  and  is  related  to  u.  The 
entropy  generated  by  the  shock  is  constant  with  respect  to  x  aft  of  the  shock.  The 
tangency  boundary  condition  is  represented  by  the  thin  airfoil  approximation  so  that 
the  equation  and  its  boundary  conditions  reduce  to  the  classic  TSD  formulation  in 
the  absence  of  vorticity.  The  tangency  boundary  condition  is  given  by 

(0y  -  =  y's(x)  (6) 

7  y  =  ±o 

where  y  =  ya(x)  denotes  the  airfoil. 

Equation  (l)  with  its  boundary  condition,  Equation  (6),  and  the  vorticity 
equation,  can  be  written  in  similarity  form  to  give 


u  =  K  +  ^y  'l 

v  =  ty-  K  J 


(3) 
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(  -2  . 
■-J- 


+  v  = 


s*L 


«2  8 


3i 


y  (7  +  i)‘ 


k  r  at  [a  -  u  Tn 


(7) 


where 


v  (x,  *o)  =  k/^3  ys  (x) 


(8) 


u  =  +  P2  iL  >  v  =  <f>-  -  fx 

y  y 


(9) 


r  i-  +  *xx  =  - u  = 


—  «3  _9 


yy 


(7  +  i) 


^  -S  [(i  -  «  +)3] 

3y 


(10) 


and  ip  =  V^2  <p;  p  =  u  =  k/p3  to,  y  =  fly  and  u  +  is  the  value  of 
u  just  ahead  of  the  shock.  In  the  derivation  of  Equation  (10)  a  first  order 
approximation  to  Crocco's  Equation  has  been  used;  thus 


u  = 


7M„ 


2  9y 


(S/R) 


(ID 


If  rys  =  yg  where  r  is  the  thickness  of  the  airfoil  then  a  transonic  similarity 


parameter,  K,  can  be  defined  as 


K  =  kT/p3  (12) 

The  transonic  similarity  parameter,  K,  is  constant  for  combinations  of  Mach 
number  and  thickness;  if  the  airfoil  is  very  thin  the  Mach  number  approaches  unity 
and  the  terms  due  to  entropy  production  vanish  as  r-»o  and  p  -*o.  Thus  for  a 
given  K  it  is  possible  to  approach  the  TSD  equation  arbitrarily  closely  by  varying 
the  thickness  and  the  Mach  number. 


A- 5 


The  set  of  equations,  Equations  (7-10),  can  be  solved  using  the  same  algorithm 
and  treatment  of  the  boundary  conditions  used  in  TSD  example,  thus  ensuring  that 
a  phantom  solution  can  be  obtained  as  and  r  approach  zero  since  in  this  case  a 
potential  solution  is  recovered. 

If  any  vorticity  is  introduced  into  the  two-dimensional  flow  the  vorticity  vector 
will  be  aligned  with  the  spanwise  direction  and  the  magnitude  of  the  vorticity  will 
not  change  in  this  direction.  According  to  the  equations  of  Klopfer  and  Nixon  (8) 
this  vorticity  is  convected  like 


(wu)x  +  (wv)y  =  0  (13) 

to  a  first  approximation.  The  velocity  induced  by  the  vorticity  is  given  by 
Equation  (4)  or  Equation  (10).  In  this  model  vorticity  is  defined  as 


v  -  u  =  u  (14) 

x  y  v  7 

Using  Equation  (14)  in  Equation  (13)  gives  an  alternative  representation  of  the 
vorticity  transport,  namely, 


xx 


A  =  u4 

yy 


u  v 
X 


u  u 

y 


(15) 


where 


A  =  w  v:  A  =  -  w  u 
x  ’  y 


(16) 
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A  suitable  boundary  condition  for  Equation  (15)  is  that  normal  derivatives  of  A  are 
specified  on  the  boundaries,  that  is,  Neumann  boundary  conditions  are  used. 


Equation  (15)  can  be  solved  using  a  simple  central  differencing  scheme  with 
the  vorticity  terms  on  the  right  hand  side  set  retarded  by  one  iteration.  This 
scheme  is  adequate  for  the  small  values  of  vorticity  introduced  in  the  present  study 
Equation  (15)  can  be  solved  iteratively  with  Equations  (l),  (3),  (4),  and  (5). 

In  some  of  the  examples  discussed  in  the  next  section  a  small  amount  of 
vorticity,  comparable  in  magnitude  to  the  shock  generated  vorticity,  is  introduced 
over  part  of  the  chord.  The  amount  of  vorticity,  u,  is  symmetric  on  both  surfaces 
of  the  airfoil. 

The  numerical  method  used  is  nearly  the  same  for  all  cases  and  is  based  on 
the  TSFOIL  [9]  code,  which  solved  the  TSD  equation  in  similarity  form.  Further 
details  can  be  found  in  Reference  [6]. 

The  TSD-Euler  equation  was  used  to  predict  the  flow  over  a  NACA  0012 
airfoil  at  M®=0.85,  a  =  0  deg,  and  at  Moo  =  0.8  and  a  =  1  deg.  The  results  are 
shown  in  Figures  (2)  and  (3),  respectively,  and  are  compared  to  solutions  of  the 
TSD  equation  and  the  Euler  equations.  It  may  be  seen  that  the  nonlifting  solution, 
Figure  (2)  compares  quite  well  with  the  Euler  solution.  The  lifting  solution  shown 
in  Figure  (3)  is  not  in  such  good  agreement.  The  probable  reason  is  that  the  thin 
airfoil  boundary  conditions  are  not  adequate.  However,  there  is  a  considerable 
improvement  over  the  TSD  equation.  It  is  interesting  to  note  that  the  dominant 
effect  of  the  entropy  is  the  change  in  shock  strength;  computations  which  included 
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only  this  effect  gave  identical  results  to  those  shown  in  Figures  (2)  and  (3).  This 
result  confirms  the  practicability  of  using  "nonisentropic"  potential  theories,  such  as 
that  by  Klopfer  and  Nixon  [10],  as  an  alternative  to  the  Euler  equations. 


Phantom  Solutions  to  the  TSD-Euler  Equation 

Before  discussing  the  occurrence  of  phantom  solution  to  the  TSD-Euler 
equation  it  is  necessary  to  indicate  how  these  solutions  were  found. 

The  form  of  the  equation  solved  is  based  on  the  TSD  formulation  solved  in 
the  code  TSFOIL;  and  is  as  follows 


12/3 

[K  -  (7+i)i  -  (7+i)iM(i>  + 1~)  +i~-K  t  - 

y  y  yy  xy 

=  ~2lL  2  M»'75  ^  <K  "  (7+1) *x  ~  C7-l)^->3  (17) 

3 (7+I) 2  x  y 


where  f  is  related  to  the  physical  perturbation  potential  $  by 


,  ,'2/3  M3/4 


=  T 


(18) 


and 
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;=r*/3Ml/2y>  K  = 


rV\ 


(19) 


The  tangency  boundary  conditions  for  $  and  ^  are 


dy 

*-  =  df  -  *-  - 0 

y  y 

where 


y  =  *  r  yg  (x) 


denotes  the  surface  of  the  airfoil. 


(20) 


(21) 


In  the  computations  a  factor  e  is  introduced  such  that  K  is  kept  constant. 
The  factor  e  is  defined  such  that 


M  t2/3  =  e  M  r  2/3  (22) 

where  Moo  and  r  are  the  new  Mach  number  and  thickness,  respectively  and  M0  and 
r0  are  values  of  the  Mach  number  and  thickness  for  which  phantom  solutions  exist 
for  the  TSD  equation.  Using  Equation  (22)  and  keeping  K  constant  gives  the 
following  relations. 


Mw  =  [  1  -  e  (1  -  M2)  ]J/2 


(23a) 


r 


r 

o 


[  £  M0/MJ3/2 


(23b) 
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In  all  of  the  computations  discussed  in  this  section  M»,  rG  remain  fixed  for  a  given 
airfoil  and  e  is  changed  to  increase  the  "Euler"  effect.  It  may  be  noted  that  if  e 
is  unity  then  M0  and  rG  are  equal  to  the  "potential"  values  M»  and  r.  A  value 
of  e  that  approaches  zero  means  that  the  TSD-Euler  equation  approaches  the  classic 
TSD  equation,  thus  ensuring  that  phantom  solutions  can  be  found  in  the  limit. 

The  grid  is  completely  symmetric  and  is  the  same  as  that  used  in  earlier 
studies  of  the  TSD  equation.  The  boundary  conditions  at  inflow  are  symmetric 
except  when  lift  is  present;  in  this  case  the  code  imposes  a  far  field  circulation  to 
speed  convergence.  This  is  identical  to  that  used  in  earlier  studies  [11]  of  the  TSD 
equation  and  does  not  influence  the  appearance  of  multiple  solutions. 

The  airfoil  geometry  for  both  sample  airfoils  is  given  in  the  code  by  an 
analytic  formula  and  symmetry  is  imposed  directly. 

Vorticity  is  introduced  through  the  surface  on  several  of  the  examples.  The 
vorticity  transport  is  modeled  by  Equation  (13)  and  the  vorticity  flux  at  the 
boundary  is  added  between  about  10%  and  45%  of  chord  according  to  the 
distribution 

«  =  {  1  +  cos  [(ir  (i-is)/|i-is|)]>  F  (24) 
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where  is  is  a  specified  point  (in  most  cases  at  x  =  0.265)  and  i  is  the  number  of 
grid  points  over  which  the  vorticity  is  introduced;  ia  is  the  mean  value  of  the 
indices  between  which  the  vorticity  is  introduced.  The  factor  F  is  introduced  to 
provide  a  skewness  to  the  distribution  in  an  attempt  to  model  the  physical  flow 
better  and  is  defined  by 


F  =  0  i  <  i 

s 

F  =  exp  (- | i-ig I )  i  > 

A  negative  value  of  w0  indicates  an  entropy  source. 


(25) 


Results 

Two  generic  airfoil  sections  were  used  in  the  study,  namely  a  NACA  00XX 
section  and  a  biconvex  section.  The  parameter  K  given  by  Equation  (19)  for  the 
NACA  00XX  section  is  kept  constant  at  the  value  for  Moo  =  0.85  and  r  =  0.12,  a 
condition  for  which  phantom  solutions  to  the  TSD  equation  can  be  found.  For  the 
biconvex  section  K  is  kept  constant  at  the  value  for  Moo  =  0.81  and  r  =  0.1.  In 
the  following  descriptions  of  the  computations  the  parameter  e  is  used  to 
differentiate  the  results;  the  corresponding  values  of  Moo  and  r  for  the  NACA  00XX 
and  biconvex  sections  can  be  found  in  Table  I  and  Table  II,  respectively.  The 
phantom  solutions  are  initiated  by  giving  the  airfoil  an  angle  of  attack  for  some 
iterations  before  putting  the  angle  of  attack  to  zero. 
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In  an  earlier  paper  it  was  noted  that  phantom  solutions  could  not  be  found 
for  this  airfoil  unless  the  Mach  number  was  close  to  unity  (actually  £  =  0.1)  but 
by  the  introduction  of  vorticity,  wQ,  the  Mach  number  range  of  the  occurrence  of 
phantom  solutions  could  be  improved.  In  the  present  study  values  of  e  between 
zero  and  0.6  were  used  and  it  was  found  that  from  a  conventional  start,  that  is, 
the  computation  is  started  as  described  earlier  from  freestream,  phantom  solutions 
could  not  be  found  for  e  >  0.2  even  when  vorticity  is  introduced;  for  e  =  0.2  the 
value  of  wQ  is  -0.125.  However,  if  the  starting  sequence  was  altered  such  that  the 
initial  flow  is  not  freestream,  then  it  proved  possible  to  obtain  phantom  solutions 
for  e  =  0.4  with  -0.105  <  uQ  <  -0.145.  In  these  results  the  lift  coefficient  was 
approximately  constant  at  a  value  of  0.33.  Examples  of  phantom  solutions  for  e  = 
0.2  and  e  =  0.4  are  shown  in  Figures  (4)  and  (5),  respectively. 

The  fact  that  the  occurrence  of  phantom  solutions  depends  on  the  starting 
solutions  indicates  that  there  may  be  more  such  solutions  than  would  commonly  be 
estimated  using  a  conventional  start  to  the  computation. 

A  phantom  solution  for  the  biconvex  airfoil  is  shown  in  Figure  (6);  the  value 
of  e  is  0.04  and  uQ  =  -0.01.  In  general  it  proved  more  difficult  to  obtain  phantom 
solutions  for  the  biconvex  airfoil. 

The  introduction  of  vorticity  through  the  boundary  has  two  effects  on  the 
computation,  namely  it  introduces  entropy  into  the  flow  field  and  it  induces  a 
velocity  on  the  airfoil  surface.  The  induced  velocity  on  the  airfoil  surface  in  a  thin 
airfoil  formulation  is  equivalent  to  changing  the  airfoil  shape.  Consequently  it  is  of 
interest  to  determine  if  the  dominant  effect  of  the  vorticity  introduction  in  the 
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appearance  of  phantom  solutions  is  the  change  in  the  effective  airfoil  geometry  or 
the  entropy  produced  in  the  flow  field.  Therefore,  the  computer  code  was  run  with 
the  change  in  geometry  being  only  effect  of  the  vorticity  introduced  through  the 
surface.  It  is  found  that  this  code  gave  essentially  the  same  result  as  the  original 
code  indicating  that  the  dominant  effect  of  the  vorticity  introduction  is  the  change 
the  airfoil  geometry  and  suggesting  that  airfoils  could  be  designed  to  give  phantom 
solutions. 

Test  for  Phantom  Solutions 

The  TSD-Euler  equation  [5]  can  be  written  in  the  form 

'*«  **„  -  -  h  •  <>&  *  *  (i  ♦  ^5)  <“) 

where 

2  ku^ 

ty  =  P‘  u  ~  2  >  K  =  -v  W 

2 

and  p  and  k  are  given  by  Equation  (2). 

If  the  substitution 

y  =  py  (28) 
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is  introduced  then  Equation  (24)  becomes 


=  ~P 

yy  3y 


(l  +  ^5)  s/R} 

tm: 


(29) 


Using  fin  analysis  similar  to  that  in  Reference  (5),  Equation  (27)  can  be  written  in 
integral  form  over  the  flow  field  D  to  give 


-  ~2  _  (1 

u  ~  I"  “  S  =  J  <A o t 


C  K 


40  s  Koy  >  d«  -  f  4o  [  1-  * 

1  °  7Moo 


♦  -4-n  a$d, 

7“. 


where 


(30) 


K  =  ^  ln)(x-02  -  (y~»7)2 1 

Ko  =  inKx-02  -  72!  . 


(31) 


(32) 


and 


Aof  =  2  (x’+°)  "  f  (x’“°)l 

Af  =  |  [f  (x , y)  -  f  (x,  -  y)] 


(33) 
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Equation  (28)  is  valid  for  y  /  0:  on  y  =  o  only  the  symmetric  part  of  Equation 
(28)  exists  and  the  asymmetric  component  is  given  by 


I 

1  —2  2 

-  ’  O'.  *o)|  -  K„x  4o  -  JJ  V  5  *  5  +  V1  ^  (34) 

average  o  '  v  '  7M00 

For  the  case  of  phantom  solutions  when  the  airfoil  is  symmetric  and  at  zero  angle 
of  attack.  V|average  is  zero  and  Equations  (28)  and  (32)  can  be  arranged  to  give 
equations  for  Au;  thus 


2  2  - 

Au  (1-u)  =  -  ff  K  -  [Auu  -  A  u  +  (1  +  -i-)  AS  -  ^  AS]  d^7 

K  VV  IK  lK 


l! 


2  _  2 

K  _  [  Auu  -  A  u  +  (1  +  -A;)  AS - ^  A  S]  dfdrj  =  0 

°  iK  iK  ° 


O  -X.TJ 


(35) 


(36) 


where  a  barred  quantity  is  defined  by 


f  =  2  [f  (x>y)  +  f  (x»_y)] 


(37) 


and  is  the  upper  half  plane. 


S  is  related  to  u  through  Equation  (5)  and  Equation  (29).  Thus  Equations 
(33)  and  (34)  are  equations  for  Au  in  terms  of  the  average  value  of  u,  u.  The 
integrals  can  be  discretized  by  a  sum;  thus  the  integral  in  Equation  (33)  can  be 
represented  by 


-lfKV ygd^Ai.gj;  i=Nx,  Nt  (38) 

while  the  integral  in  Equation  (34)  can  be  represented  by 

if  KX,  **«<">  =  ®ij  *i  1  -  »-  ",  (39> 

l 

where  Nx  and  Ny  and  the  number  of  grid  points  in  the  x  and  y  directions 
respectively  and 


Nt  =  NxNy  (40) 

Since  AS  is  a  function  of  Au  and  u  Equations  (33)  and  (34)  can  be  represented  by 
the  following  matrix  equation. 


Aij  =  0  i  =  1  Nt  (41) 

Note  that  both  Equation  (33)  and  Equation  (34)  have  been  combined  to  give  this 
result.  Note  also  that  A;j  is  a  function  of  u  and  Au.  For  the  potential  TSD 
formulation  A;j  is  a  function  only  of  u. 

Suppose  that  a  phantom  solution  can  occur  with  an  infinitesimal  amount  of 
asymmetry  that  is  Au  can  be  vanishingly  small;  evidence  for  Au  being  vanishingly 
small  can  be  inferred  from  the  results  of  Williams  et  al(4)  where  the  phantom 
solution  can  have  infinitesimal  lift.  In  such  a  case  the  symmetric  solution  bifurcates 
at  some  point  and  at  this  bifurcation  point  both  symmetric  and  asymmetric 
solutions  are  possible.  This  can  only  occur  if  the 
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det  A. .  =  0 
ij 


(42) 


where  "det  ()"  denotes  the  determinant  of  the  matrix  Ajj  at  this  point.  Since  at 
this  point  Au  is  vanishingly  small  A;j  is  a  function  of  u  and  hence  Equation  (40)  is 
an  equation  for  u.  Since  in  the  TSD  approximation  u  is  related  to  the  local  Mach 
number  distribution  (M)  by 


u 


M 


m: 


(43) 


Equation  (40)  is  an  equation  for  the  Mach  number  distribution  in  terms  of  the 
freestream  Mach  number.  Hence,  the  occurrence  of  phantom  solutions  is  dependant 
only  on  the  Mach  number  distribution  for  a  fixed  M».  In  other  words  there  is  a 
Mach  number  distribution  for  which  the  solution  can  bifurcate,  indicating  that 
phantom  solutions  can  be  "designed"  to  appear.  This  fact  is  borne  out  by  the 
results  of  the  vorticity  introduction  study  which  indicated  that  a  dominant  cause  of 
phantom  solutions  is  the  airfoil  geometry.  Thus,  it  should  be  possible  to  design  an 
airfoil  which  will  have  phantom  solutions  for  a  given  Moo. 


In  order  to  test  whether  the  phantom  solution  satisfy  Equation  (36)  the 
phantom  solution  for  the  Joukowski  airfoil  in  potential  flow  at  Moo  =  0.85  can  be 
used  to  evaluate  the  integral.  However,  evaluation  of  the  integral  is  not 
satisfactory,  principally  due  to  the  fact  that  the  shock  is  captured  and  thus  is  not 
a  solution  of  the  differential  equation,  furthermore,  at  the  shock  the  algorithm  has  a 
truncation  error  of  0(1).  A  better  means  of  judging  whether  the  integral  is  zero  is 
to  used  an  inverted  form  of  Equation  (36),  namely, 
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1  r  Au2  1  1  1  -  X  1/2  f1  (  ^2 

F(x)  =  \  [  Au(x)  -  ^2~  (x)]  +  \  J  7T7  (r^) 


d£  =  0  (44) 


where,  for  the  TSD  equation, 


^0  =  U  ^oxt]  [  Auu  -  Aouuo]  d£d7 
U1 


(45) 


In  this  formulation  the  integral  Ic(£)  is  not  so  sensitive  to  errors  since  on  tj 
=  0  the  integrand  is  zero  and  the  increasing  error  as  r/  «  is  counteracted  by  the 
decreasing  of  the  kernel  function  Kox^  as  y+  Since  for  a  normal  shock  the  first 
term  is  a  smooth  function  the  error  in  the  shock  capture  region  can  be  reduced  by 
interpolation  of  this  term  for  points  inside  the  shock  capture  region.  As  noted 
earlier,  the  integral  Ic,  and  its  inverse,  should  not  be  too  inaccurate.  The  variation 
of  F(x)  and  x  is  shown  in  figure  7.  The  value  for  the  Joukowski  airfoil  at  Mo#  = 
0.81  and  a  =  1/4*  and  its  exact,  theoretical  solution,  are  shown  for  comparison.  It 
can  be  seen  that  the  F(x)  is  very  close  to  zero  (apart  from  the  shock  region)  and 
strongly  suggests  that  the  phantom  solutions  satisfy  Equation  (36).  Evaluations  of 
the  integral  for  other  configurations  lead  to  the  same  conclusion.  Thus,  while  it  is 
unlikely  that  the  finite  difference  results  will  give  the  singular  matrix,  det  A;j,  at 
the  same  freestream  conditions  as  the  actual  occurrence  in  the  calculations.  The 
determinant  should  become  singular  at  approximately  the  correct  Mach  number. 
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In  order  to  test  the  criterion  of  Equation  (40)  further  the  determinant  of  Ay 
was  evaluated  by  performing  a  LU  decomposition  of  the  matrix  and  identifying  the 
smallest  magnitude  of  the  elements  on  the  diagonal.  For  a  singular  matrix  this 
magnitude  should  be  zero.  As  a  test  a  true  singular  equation,  based  on  Equations 
(33)  and  (34)  was  set  up  by  multiplying  the  integrals  be  a  factor  such  that  these 
equations  are  satisfied  identically  by  a  finite  difference  phantom  solution.  The 
smallest  element  on  the  diagonal  of  the  LU  decomposition  is  then  taken  to  be  zero 
in  subsequent  cases.  Note  that  the  smallest  element  is  not  zero  because  of  machine 
error  and  the  overall  magnitude  of  the  determinant;  in  the  present  study  the 
determinant  itself  could  be  as  small  as  exp (-200). 

For  the  TSD  equation  the  calibration  test  noted  earlier  gives  a  minimum 
magnitude  of  5.64  x  10"4  and  in  the  subsequent  tests  anything  less  than  this  value 
is  assumed  to  be  (probably)  zero.  The  machine  used  is  a  VAX  11/750  computer. 
In  Figure  (8)  the  magnitude  of  the  minimum  diagonal  element  is  shown  for  various 
values  of  Moo  for  NACA  0012,  14%  biconvex  and  11%  Joukowski  airfoil.  It  can  be 
seen  that  phantom  solutions  can  be  found  in  the  neighborhood  of  Moo  “  0.81  for 
the  NACA  0012,  Moo  ~  0.88  for  the  Joukowski  airfoil  and  Moo  ~  0.86  for  the 
biconvex  airfoil,  Moo  ~  0.86.  The  range  for  multiple  solutions  for  the  NACA  0012 
and  11.8%  Joukowski  airfoils  using  the  TSD  equations  is  0.82  -  0.86,  and  for  the 
14%  biconvex  airfoil,  0.79-0.83. 

Although  the  test  appears  to  work  there  are  some  aspects  of  the  derivation  of 
the  test  that  deserve  a  warning.  First,  the  basis  of  the  test  is  that  there  is  a 
bifurcation  point  where  the  symmetric  solution  can  change  smoothly  to  a  phantom 
solution.  It  is  not  proven  that  such  a  bifurcation  point  exists,  although  the  results 


of  Reference  [4]  do  indicate  such  a  bifurcation.  Secondly,  the  inability  of  the 
present  formulation  to  give  machine  zero  (10’6)  for  a  real  singular  equation  can  cast 
doubt  on  the  conclusions. 

The  above  tests  are  for  the  TSD  equation.  A  similar  series  of  tests  were 

computed  for  the  TSD-Euler  equation.  In  these  tests  the  parameter  e  in  Equation 

(23)  is  picked  so  that  the  airfoil  thickness  is  constant  for  a  range  of  Mach 

numbers.  The  thickness  is  picked  so  that  at  Mx  =  0.986027  there  is  a  phantom 

solution.  Results  of  the  test  are  shown  in  Figure  (9).  The  "zero"  obtained  in  a 

similar  fashion  as  for  the  TSD  equation  for  a  grid  similar  to  that  used  in  the  TSD- 

-3 

Euler  calculations  is  2.1  x  10’  .  It  may  be  seen  from  Figure  (9)  that  the  test  does 
not  work  as  well  for  the  TSD-Euler  equation  as  for  the  TSD-Equation  although  it 
is  possible  that  there  is  a  bifurcation  at  M®  =  0.9753  which  gives  a  value  close  to 
the  "zero"  value. 

The  test  described  above  is  logical  in  a  mathematical  sense  and  can  give 
insight  into  the  occurrence  of  phantom  solutions.  However,  as  a  predictor  or  a 
design  criteria  for  phantom  solutions  it  requires  further  work  to  remove  the 
obstacles  of  numerical  error  accumulation  and  the  fact  that  the  numerical  solutions 
do  not  satisfy  the  integral  equations  that  are  the  basis  of  the  theory. 
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Multiple  Solutions  of  the  Navier-Stokes  Equations 

The  type  of  phantom  solutions  discussed  in  the  previous  sections  have  been 
solutions  of  the  TSD  or  TSD-Euler  equation.  However,  it  is  known  that  multiple 
solutions  of  the  Navier-Stokes  equations  can  occur  at  transonic  speeds.  For  example 
the  results  of  Levy  [12]  showed  that  a  18%  biconvex  airfoil  at  zero  angle  of  attack 
will  give  an  oscillatory  solution  for  a  narrow  range  of  Mach  numbers.  Other 
solutions  of  this  type  have  been  discussed  by  Mabey(l3)  and  a  range  of  thickness 
parameters  for  which  these  multiple  solutions  appear  has  been  determined.  The 
occurrence  of  these  solutions  are  not  limited  to  the  biconvex  airfoil.  The  term 
multiple  solutions  is  used  here  because  both  the  symmetric  and  the  oscillatory 
solutions  have  been  observed  in  wind  tunnel  tests  so  the  terminology  "phantom" 
solution  may  be  misleading.  In  the  present  study  the  computer  code  ARC2D  [14] 
was  used  to  computer  these  oscillatory  solutions  for  a  14%  biconvex  airfoil  at  zero 
angle  of  attack  and  M»  =  0.81;  the  Reynolds  number  is  6.9  x  10  .  The  trace  of 
the  lift  coefficient  with  time  is  shown  in  Figure  (10). 

The  range  of  Mach  numbers  in  which  these  oscillatory  solutions  appear  is  close 
to  that  for  which  phantom  solutions  appear  in  the  TSD  equation  and  it  is  a 
natural  question  as  to  whether  the  oscillatory  solution  for  the  Navier-Stokes 
equations  in  initiated  by  a  mechanism  similar  to  that  for  the  inviscid  flows.  The 
oscillation  is  probably  sustained  by  viscous  interaction.  Accordingly  a  series  of 
calculations  were  performed  to  try  to  stabilize  the  asymmetric  solution.  This  was 
attempted  by  a  symmetric  boundary  layer  suction  and  by  attempting  a  solution  for 
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a  very  thin  airfoil  at  M»  close  to  unity  to  mimic  the  results  for  the  TSD-Euler 
equation.  The  suction  failed  to  produce  a  stable  solution  while  the  very  thin  airfoil 
flow  became  dominated  by  the  boundary  layer  and  attempts  to  remove  this  by 
suction  were  not  successful. 

Although  the  attempts  to  produce  a  stable  solution  for  the  Navier-Stokes 
equations  failed  it  is  necessary  to  note  that  the  Navier-Stokes  equations  do  have 
multiple  solutions  in  a  small  Mach  number  range  since  there  is  one  symmetric  and 
one  oscillatory  solution. 

If  the  Navier-Stokes  equations  are  regarded  in  the  same  light  as  an  inviscid 
flow  with  a  boundary  layer  it  is  possible  to  construct  a  scenario  for  the  appearance 
of  the  multiple  solutions.  If  the  boundary  layer  can  be  represented  by  a 
displacement  thickness  then  the  symmetrical  flow  at  the  start  of  the  oscillatory  cycle 
can  be  regarded  as  an  inviscid  flow  over  a  modified  airfoil  section.  Alternatively 
the  flow  outside  a  boundary  which  borders  the  boundary  layer  can  be  regarded  as 
inviscid.  Since  it  has  been  shown  that  the  effect  of  entropy  as  regards  the 
appearance  of  phantom  solutions  is  mainly  to  modify  the  Mach  number  distribution 
this  inviscid,  nonisentropic  flow  can  give  phantom  solutions.  Suppose  now  that  at 
the  start  of  the  oscillatory  cycle  a  phantom  solution,  in  the  inviscid  sense,  occurs. 
This  solution  will  create  circulation  and  drive  one  shock  wave  aft  and  the  other 
forward.  Since  the  shock  movement  will  be  associated  with  a  change  in  shock 
strength  the  boundary  layer  separation  region  will  change  in  both  extent  and 
strength  causing  a  change  in  displacement  thickness  which  changes  the  effective 
shape  of  the  airfoil.  This  modified  airfoil  could  be  outside  the  range  of  phantom 
solutions  so  the  solution  will  then  tend  towards  symmetry.  Since  this  is  a  dynamic 
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process  an  oscillatory  cycle  could  be  produced.  This  scenario  is  really  only 
speculation  but  some  evidence  for  its  existence  can  be  provided  by  a  simple  time 
dependent  computation. 

Consider  the  potential  flow  over  a  symmetric  airfoil  at  zero  angle  of  attack  in 
the  range  of  phantom  solutions.  A  small  perturbation  to  the  boundary  is  given  so 
that  a  phantom  solution  starts  to  be  generated.  As  this  happens  one  shock  will 
get  stronger  and  a  viscous  flow  would  separate  the  boundary  layer.  The  added 
displacement  thickness  from  this  separation  can  be  represented,  at  least 
approximately,  by  the  "viscous  wedge"  model  suggested  by  Rizzetta  &  Yoshihara 
[15].  In  the  present  formulation  the  airfoil  slope  is  changed  by  the  addition  of  the 
term  af  *  where 
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where  xs  is  the  shock  location  and  8  represents  the  slope  of  the  separation 
displacement  thickness,  given  by 


8  = 


0,  M  <  1.3 

8  (M+  -  1),  M+  >  1.3 


(47) 


M+  is  the  Mach  number  just  ahead  of  the  shock;  8Q  is  to  be  determined. 
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The  computer  code  XTRAN2L  [14]  was  modified  to  incorporate  the  viscous 
wedge;  XTRAN2L  solves  the  time  accurate  TSD  equation.  The  airfoil  used  in  the 
test  was  a  14%  biconvex  airfoil  at  Moo  =  0.81,  the  same  as  in  the  Navier-Stokes 
computations.  The  phantom  solution  was  obtained  by  an  initial  perturbation  and 
the  shock  strength  on  the  upper  surface  was  such  that  M+  >  1.3  thus  triggering 
the  viscous  wedge.  Several  values  of  6Q  were  tried  and  the  response  of  to  the 
shock  was  lagged  in  time.  In  order  to  reproduce  qualitatively  the  Navier-Stokes 
solution,  0O  was  taken  to  be  0.32,  considerably  larger  than  that  recommended  by 
Rizzetta  and  Yoshihara  [14]  (0O  =  0.09);  the  viscous  wedge  was  introduced  50  time 
steps  behind  the  shock  motion.  In  the  subsequent  oscillation  this  is  equivalent  to  a 
lag  of  approximately  a  quarter  of  a  cycle.  The  oscillation  of  the  lift  coefficient  is 
shown  in  Figure  (11).  Only  part  of  the  cycle  is  shown;  the  cycle  persists  for  the 
complete  calculations,  consisting  of  1500  time  steps.  It  can  be  inferred  from  these 
computations  that  a  plausible  initial  cause  of  the  oscillatory  solutions  of  the  Navier- 
Stokes  equations  is  the  same  as  that  for  the  inviscid  cases  but  that  viscous  forces 
are  responsible  for  the  oscillatory  behavior.  It  is  important  to  remember  that  the 
Navier-Stokes  solutions  do  occur  in  wind  tunnel  tests  and  are  not  a  mathematical 
or  computational  anomaly. 


Concluding  Remarks 

The  TSD-Euler  equation,  derived  in  previous  work  was  examined  for  phantom 
solutions  in  an  attempt  to  determine  the  effect  on  entropy  and  vorticity  on  the 
appearance  of  such  solutions.  It  was  found  that  phantom  solutions  can  be  found 
for  this  equation  although  for  a  different  range  of  Mach  numbers  than  for  the  TSD 
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equation,  at  least  for  the  airfoils  considered.  Following  earlier  work  a  source  of 
vorticity  was  introduced  on  the  airfoil  surface  and  this  enhanced  the  probability  of 
the  appearance  of  phantom  solutions.  Closer  study  indicated  that  the  dominant 
effect  of  the  vorticity  introduction  is  to  change  the  effective  geometry  of  the  airfoil 
with  the  implication  that  an  airfoil  could  be  designed  to  give  a  phantom  solution  to 
the  TSD-Euler  equation.  The  vorticity  introduction  can  be  regarded  as  a  method 
of  flow  control  which  mimics  the  effect  of  a  change  in  geometry. 

A  test  of  when  a  phantom  solution  can  appear  has  been  developed  and  it 
indicates  that  the  dominant  feature  determining  the  appearance  of  phantom  solutions 
is  the  local  Mach  number  distribution;  the  entropy  or  vorticity  in  the  flow  field 
plays  only  a  secondary  role.  The  numerical  verification  of  the  test  has  some 
problems  due  principally  to  the  numerical  inaccuracies  in  the  numerical  solutions  of 
the  flow  and  in  the  evaluation  of  the  test  parameters.  The  test  indicates  that  a 
phantom  solution  can  be  initiated  by  designing  the  airfoil  to  give  a  suitable 
distribution  of  Mach  number  in  the  flow  Field. 

The  inference  of  the  study  of  the  TSD-Euler  equation  is  that  entropy  and 
vorticity  do  not  inhibit  the  appearance  of  phantom  solutions  but  simply  change  the 
freestream  conditions  for  which  these  solutions  can  appear.  Since  vorticity  does  not 
inhibit  the  appearance  of  the  phantom  solutions  it  is  pertinent  to  ask  if  the  Navier- 
Stokes  equations  have  phantom  solutions.  It  is  known  that  for  a  small  Mach 
number  range  that  the  Navier-Stokes  equations  do  have  multiple  solutions,  are 
symmetric  and  are  oscillatory.  This  behavior  was  shown  to  be  explainable  by  the 
idea  of  phantom  solutions,  although  this  may  not  be  the  only  explanation.  The 
multiple  solutions  of  the  Navier-Stokes  equations  are  apparent  in  wind  tunnel  tests, 
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indicating  their  realism.  The  question  that  now  arises  is  whether  phantom  solutions 
can  be  made  to  exist  stably  in  a  real  environment  thus  introducing  a  new  area  of 
study  for  aerodynamicists  and,  possibly,  introducing  a  new  mechanism  for  flight.  In 
order  to  achieve  this  the  phantom  solutions  must  be  stable  and  some  form  of  flow 
control  is  necessary. 

If  the  test  derived  in  this  present  work  can  be  made  more  accurate 
numerically  it  should  be  possible  to  design  an  airfoil  to  give  phantom  solutions  over 
a  range  of  flight  conditions.  Again,  some  form  of  flow  control  is  necessary.  If 
such  an  airfoil  can  be  designed  then  phantom  solutions  to  the  Navier-Stokes 
equations  may  be  found,  together  with  some  means  of  stabilizing  such  solutions. 
After  this  stage  a  wind  tunnel  test  could  be  conducted.  To  achieve  such  an  end  it 
is  necessary  to  learn  much  more  about  the  nature  of  all  phantom  solutions,  not  just 
those  considered  here,  that  is,  the  transonic  solutions. 
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Figure  8a.-  Test  for  p!„-nuoi  solutions  lor  TSD 
equation;  NAt.'A  0012  airfoil. 
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Figure  9.-  Tesl  for  phantom  solutions  for  TSD- 
Euler  equation;  NACA  OOXX  airfoil. 


A-37 


0.0  300.0  600.0  900.0  IT 

ITER 

Figure  JO.-  Trace  of  vs.  time  tor  14/.',  biconvex 
airfoil,  M  =  0.81,  a  =  0.0  deg;  Navier- 


Stokes  solution. 
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